<p> Draghilev RK </p>

8c112f66-3f68-4539-bb4c-ab4f7451636c 10 6 9 decimal

&[DATE] &[TIME] - &[FILENAME]

X0 tmin tmax N D n f rows 1 + : f n el 0 : j 1 : j n ≤ j j 1 + : k 1 : k n ≤ k k 1 + : m j k el f j el x k el diff : for for Матрица Якоби x Jacobi m : Матрица  правых частей алгебраических уравнений b 1 - (x Jacobi 1 n 1 - n n submatrix * : Определители A x Jacobi 1 n 1 - 1 n 1 - submatrix : n 2 > Δ 1 el b A 1 n 1 - 2 n 1 - submatrix augment det - : Δ n 1 - el A 1 n 1 - 1 n 2 - submatrix b augment det - : k 2 : k n 2 - ≤ k k 1 + : Δ k el A 1 n 1 - 1 k 1 - submatrix b A 1 n 1 - k 1 + n 1 - submatrix augment det - : for Δ n el A det - : 41 line Δ 1 el b : (Δ n el A det : (& if i 1 : i n ≤ i i 1 + : D2 i el Δ i el Δ s el (2^(s 1 n sum sqrt / : for t x D1 D2 : result X0 tmin tmax N t x D1 rkfixed 1 N 2 X0 rows 1 + submatrix : 131 line :

2.Continuously Stirred Tank Reactor

http://demonstrations.wolfram.com/ContinuouslyStirredTankReactorUsingArcLengthParameter/

x1-dimensionless concentration, x2-dimensionless temperature, Da-Damköhler number, γ-dimensionless activation energy, β-dimensionless heat transfer coefficient

B 10 :

Da 0.03084 :

γ 20 :

T0 5 - :

Kоординаты начальных точек

Denote that Smath 5346 supports the "compactum style"

tmin 0 : tmax 9.8 - : Δt 0.2 - : N tmax Δt / : 41 line

X0 0.004 1.85 - 0.03084 31 mat :

x x T Da 31 mat ≡

The dimensionless equations

φ x 2 el - B x 3 el * 1 x 1 el - (* e x 2 el 1 x 2 el γ / + /^* + :

f 1 el x 1 el - x 3 el 1 x 1 el - (* e x 2 el 1 x 2 el γ / + /^* + : 0 ≡

tstart 6 :

tend 30 :

i tstart tend range β i el 0.2 0.04 i * + : f 2 el φ β i el x 2 el * - : K i el X0 tmin tmax N D : 31 line for

case insensitive ... very poor Greek !
Takes only Greek "tau"  [ι] ι tstart tend range :

Dependence T of on the parameter Da. By increasing the heat transfer coefficient the temperature drops

t tend ≥ t 2 tend * t - : t t : if K t el 3 col K t el 2 col augment 2 1 sys

tt 6 :

Export K tt el 3 col K tt el 2 col augment :

tt from 6 ... 30

Export

multisol K 6 el 3 col K 6 el 2 col augment K 30 el 3 col K 30 el 2 col augment K 12 el 3 col K 12 el 2 col augment K 18 el 3 col K 18 el 2 col augment K 24 el 3 col K 24 el 2 col augment 51 sys :

multisol